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ABSTRACT 

Glint, the specular reflection of sunlight off Earth's oceans, may reveal the 
presence of oceans on an extrasolar planet. As an Earth-like planet nears crescent 
phases, the size of the ocean glint spot increases relative to the fraction of illumi- 
nated disk, while the reflectivity of this spot increases. Both effects change the 
planet's visible reflectivity as a function of phase. However, strong forward scat- 
tering of radiation by clouds can also produce increases in a planet's reflectivity 
as it approaches crescent phases, and surface glint can be obscured by Rayleigh 
scattering and atmospheric absorption. Here we explore the detectability of glint 
in the presence of an atmosphere and realistic phase- dependent scattering from 
oceans and clouds. We use the NASA Astrobiology Institute's Virtual Planetary 
Laboratory 3-D line-by-line, multiple-scattering spectral Earth model to simulate 
Earth's broadband visible brightness and reflectivity over an orbit. Our validated 
simulations successfully reproduce phase-dependent Earthshine observations. We 
find that the glinting Earth can be as much as 100% brighter at crescent phases 
than simulations that do not include glint, and that the effect is dependent on 
both orbital inclination and wavelength, where the latter dependence is caused 
by Rayleigh scattering limiting sensitivity to the surface. We show that this 
phenomenon may be observable using the James Webb Space Telescope (JWST) 
paired with an external occulter. 

Subject headings: astrobiology — Earth — planets and satellites: composition 
— radiative transfer — techniques: photometric — scattering 
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Introduction 



A major goal in the study of extrasolar planets is the detection and recognition of 
a "habitable" world, or a planet capable of maintaining liquid water on its surface. A 
variety of direct and indirect approaches could be used to determine if a planet is habitable. 
Indirect approaches focus on characterizing the surface environment of a planet, which would 
constrain the likelihood that the planet could maintain liquid water on its surface. Direct 
appr oaches aim to detec t signs that indicate the presence of water on the surface of a planet 



(e.g. 



Cowan et al.ll2009f ). 



One direct indicator of surface bodies of water is specular r eflection, or the "glin t effect 



While specular reflection is not unique to liquid surfaces (e.g., iDumont et al.ll2010l ). liquids 
are distinguished from other surfaces by their contrast between weak specular reflectance 
at direct illumi nation angles and str ong specular reflectance at glancing illumination an- 
gles. Recently, IStephan et al.l ( 120101 ) used the enhanced reflectivity of liquids at glancing 
illumination angles to provide evidence for liquid hydrocarbon lakes on the surface of Titan. 



Sagan et al.l ( 119931 ) argued that the presence of a specularly-reflecting region (or "glint 



spot") in spatially-resolved images of Earth taken by the Galileo spacecraft, combined with 
detections of atmospheric water vapor and surface temperatures near the melting point of 
water, was evidence for the presence of liquid water oceans on Earth's surface. Unfortunately, 
obtaining spatially-resolved images of terrestrial extrasolar planets presents an engineering 
challenge that will not be met in the near future. The first measurements that aim to detect 
glint must rely on how it affects the brightness of a planet in a disk-integrated sense. 

The relative size of Earth's glint spot compared to the illuminated portion of the disk 
increases at crescent phases, and the reflectivity of water increases at glancing illumination 
an gles, affecting Earth ' s pha se curve. The detectability of this effect was first investigated 
by lWilliams fc Gaidod (120081 ) . who used a simple model of Earth's reflectance coupled to a 
3-D climate model to predict Earth's appearance over the course of an orbit. Their model 
showed that Earth's reflectivity increases into c rescent phases, b ut was unable to reproduce 
Earthshine observations of Earth's reflectivity ( IPalle et al.ll2003l ). This discrepancy was at- 
tributed to the absence of Rayleigh scattering in their model and the assumption that clouds 
reflect isotropically (Lambertian). In reality, liquid droplets and ice crystals preferentially 
scatter light in the forward direction, which can mimic the glint effect. 



Oakley fc Cash! ( 120091 ) modeled Earth's brightness over the course of an orbit with a n 



emphasis on characterization by the New Worlds Observer mission concept (ICashl 120061 ). 
Their model simulated Earth's reflectivity using satellite-measured bidirectional r eflectance 



distribution functions (BRDFs) for a variety of scenes (e.g., thick cloud over ocean) (IManalo-Smith et al 
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19981 ). and used satellite observations to evolve clouds, snow, and ice in their simulations. 
Their model also demonstrated an increase in Earth's reflectivity into crescent phases, but 
was not compared to the phase- dependent Earthshine measurements. The authors proposed 
that the bright glint spot increases the variability in Earth's brightness at crescent phases, 
which could serve as an indicator of surface oceans. However, the BRDFs used in their model 
were not valid at extreme crescent phases, requiring assumptions about cloud scattering and 
ocean reflectivity at glancing illumination angles. 

In this work we use the NASA Astrobiology Institute's Virtual Planetary Laboratory 
(NAI-VPL) 3-D spectral Earth model to simulate Earth's disk-integrated spectrum as it 
would appear to a distant observer watchi ng the planet through an orbit. While extensively 



validated against time-dependent data in iRobinson et al.l (120101 ). we present further valida- 



tions here against phase-dependent Earthshine observations. The validated model is used 
to investigate the significance of glint in Earth's phase curve by discriminating between the 
competing effects of cloud scattering and glint. We also discuss the observing requirements 
for glint detection. 



2. Model Description 



The NAI-VPL 3-D spectral Earth model simulates Earth's appearance to a distant 
observer. The model produces spatially- and spectrally-resolved datacubes that can be col- 
la psed to produce high-r esolution, disk-integrated spectra, and was described and validated 



m 



mation can be found in 



Robinson et al.l ( 12010), so only a b r ief sum mary will be presented here. Additional infor- 



Tinetti et al.l (l2006al Jbl). 



The 3-D spectral Earth model simulates Earth's spectrum at arbitrary viewing geome- 
try over wavelengths from the far-ultraviolet to the far-infrared on timescales from minutes 
to years. Earth-observing satellites provide spatially-resolved, date-specific inputs of key 
surface and atmospheric properties. To simulate the seasonal changes in Earth's appear- 
ance over a year, we use snow cover and sea ice data as well as cloud cover and optical 
thickne ss data from the Moder ate Resolution Imaging Spectroradiometer (MODIS) instru 



ments (jSalomonson et al.l 119891 ) aboard NASA's Terra and Aqua satellites (IHall et al.l 11995 
Riggs et al.lll999[ ). 



Specular reflec tance from liquid wa ter surfaces in our model is simulated using the Cox- 
Munk glint model ( ICox fc Munk! 119541 ). which allows for the calculation of the BRDF of a 
wave-covered ocean given wind speed and direction, which are provided by the QuikSCAT 
satellite (http://winds.jpl.nasa.gov/missions/quikscat/index.cfm). Wavelength-dependent 
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optical properties for liquid clouds are derived usin g a Mie theory model ( ICrispl 119971 ) and 
ice clouds are parametrized using geometric optics (jMuinonen et al.lll989l ). Figure [1] shows 
a true-color image generated by our model, demonstrating the glint effect and the ability of 
clouds to mimic this effect. 



Model Validation 



Our model was p reviously validated ag ainst time-dependent ultraviolet through mid- 
infrared observations ( IRobinson et al.ll2010l ). but for a narrow range of phases. Here we 
extended our validation by matching Earthshine m easurements of E arth's apparent albedo 
at a variety of phases between gibbous and crescent (IPalle et al.ll2003l ). The apparent albedo 
is defined as the albedo of a perfect Lambert sphere that would give the same reflectivity 
of a body at a given phase angle. Thus, the apparent albedo of a Lambert sphere would be 
constant through all phases. 

Figure [2] shows our simulation of Earth's brightness and corresponding apparent albedo 
as it would appear over the course of o ne year to a distant ob server (black curves). Obser- 
vations from NASA's EPOXI mission (ILivengood et al.l 120081 ) as well as a large number of 
ground-based Earthshine observations are also shown. Observations and model results span 
the wavelength range 0.4-0.7 /im. EPOXI observations were recorded over three separate 
24-hour periods in March, May, and June of 2008, while Earthshine measurements span 
November 1998 through January 2005. Calib ration errors as large as 10% are present in 
the EPOXI obse rvations ([Klaasen et al.l 120081 ) while Earthshine observations have a stated 
accuracy of 2% jQiu et al.ll2003h . 



In general, there is good agreement between the data and the model, demonstrating the 
ability of the model to properly simulate Earth's phase-dependent brightness. At gibbous 
phases where Earthshine data were recorded, the mean apparent albedo of the model Earth 
is 0.29 ±0.02 and is 0.26 ±0.03 for the Earthshine data. Thus, while the Earthshine data are 
systematically lower than the EPOXI observations and the model at gibbous phases, they 
are all in agreement to within one standard deviation. At crescent phases the Earthshine 
observations are systematically larger than the model. However, an analysis where the 
Earthshine data were divided into 10°-wide bins in orbital longitude shows that the model 
is always within a single standard deviation of the observations. 
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4. Results 

In the subpanel of Fig. [2J the glinting model demonstrates an excess in brightness as 
large as 50% over the non-glinting model (grey curves), which uses an isotropically-scattering 
ocean reflectance model whose albedo reproduces Earth's geometric albedo to within a few 
percent. The excess brightness increases into crescent phases as the contribution of the glint 
spot to Earth's disk-integrated brightness grows, peaking at orbital longitudes near 20-30° 
from new phase. The brightness excess declines as the ill uminated cres c ent sh rinks further. 



This is due to a wave-surface "hiding" effect discussed in lCox fc Munkl (11954 ). where ocean 
waves block rays of light at glancing illumination angles. 

The brightness difference between the glinting and non-glinting model is a strong func- 
tion of wavelength. If we instead choose a filter that spans 1.0-1.1 yum, the peak brightness 
excess increases to about 100%. This behaviour arises because glint occurs at the surface 
and some wavelength ranges are more sensitive to Earth's surface than others. 

Near 90° orbital longitude, both models in Fig. [2] have apparent albedos that are 
ab out 15-20% larger than th ose near 270°. This is due to Earth's seasons, as was noted 



in 



Williams fc Gaidod (120081 ). and was determined by comparing to a model run without 
any seasonal evolution of snow and ice. As Earth moves from northern winter to northern 
summer (0° to 180° orbital longitude in these simulations, respectively), the illumination of 
the northern polar region, which is tilted towards the observer, increases. Since this region 
of Earth is more snow- and ice-covered prior to northern summer than following northern 
summer, the planet appears more reflective heading into full phase than moving out of full 
phase. The ma g nitud e of this asymmetry in the lightcurve agrees with the simulations of 



Oakley fc Cash! (120091 ) . In general, the effects of seasons on Earth's lightcurve depend on 
viewing geometry, but are small compared to the effects of cloud scattering and glint at 
crescent phases. 

The variability of both models in Fig. |2J defined as the ratio between the standard 
deviation of all model observations from a 24-hour period and the 24-hour average brightness 
from the same timespan, is shown in Fig. [3j The variability steadily increases from about 5% 
near full phase to about 30-40% near crescent phase. Brightness is more variable at crescent 
phases since the illuminated portion of the disk represents a relatively small fraction of the 
planet's surface area and, thus, is easily dominated by clouds that rotate into view or a 
cloud- free view of the glint spot. Variability at gibbous phases following northern summer 
is slightly larger than variability at gibbous phases prior to northern summer (4% versus 
6%, respectively), which is a season al effect. The m a gnitu de of the variability in our models 



agrees well with the simulations of lOakley fc Cash! (120091 ). However, these authors did not 



find a seasonal dependence in variability measurements. Furthermore, their simulations show 
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a sharp increase in variability as the planet moves into crescent phases (variability increases 
from 5-10% to 40% over about 10° of orbital longitude) while our models show a gradual 
increase in variability into crescent phases. 



Discussion and Conclusions 



5.1. Earth With and Without Glint 



Including phase-dependent reflection from oceans and clouds as well as Rayleigh scat- 
tering has allo wed us to reproduce both the brightness and phase dependence of Earthshine 
observations. IWilliams &: Gaidosl (120081 ) explicitly ignored Rayleigh scattering and phase- 
dependent scattering from clouds and were unable to reproduce Earthshine observations, 
demonstrating the importance of including these effects in a realistic spectral Earth model. 
A model that only includes phase-dependent Rayleigh scattering produces an increase in ap- 
parent albedo at crescent phases (due to weak forward and backward scattering lobes in the 
Rayleigh scattering phase function), but the upturn occurs only at extreme crescent phases 
(at orbital longitudes within 30° of new phase), which is not seen in Earthshine observa- 
ti ons. This argues tha t the l ack of phase- dependent cloud scattering in the model presented 



m IWilliams &: Gaidosl ( 120081 ) was the primary reason why their model could not reproduce 
Earthshine observations. Thus, predictions regarding the behavior of Earth's brightness at 
crescent phases are especially reliant on realistic cloud modeling. 

Models that do not include the "hiding" effect of ocean waves will over-estimate the 
brightness of water surfaces at glancing reflection angles. For edge-on orbits (inclination, 
i = 90°), this effect becomes esp ecially important at orbital longitudes within about 30° 
of new phase, in agreement with ICox fc Munki ( 119541 ). The phase-dependent relative size 
of the glint spot, the tendency of water to be more reflective at glancing reflection angles, 
and the "hiding" effect all combine to produce a maximum brightness excess for a realistic 
Earth over an Earth without glint near 30° from new phase (for an orbit viewed edge-on). 
Varying ocean wind speeds in our model show that the location of this peak is only weakly 
dependent on s urface wind conditions. Previous models used t o investigate the detection of 
surface oceans ( IWilliams fc Gaidosl 120081 ; lOakley fc Cash! 120091 ) do not include the "hiding" 
effect and cannot make strong statements about glint detection at extreme crescent phases. 

The season- and phase-dependent variability of Earth's brightness, shown in Fig. [3] 
and taken from the edge-on simulations shown in Fig, [2J is due to contrast between highly 
reflective surfaces and surfaces with low reflectivity. Following full phase, which corresponds 
to northern summer in our simulations, snow and sea ice in the northern polar region have 
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been replaced by darker surfaces (e.g., grassland) which provide greater contrast to clouds as 
compared to the snow and ice present prior to full phase. Thus, variability is larger following 
northern summer in our simulations. At crescent phases, contrast is provided by bright, 
forward-scattering clouds, and/or the bright glint spot, against relatively dark, Lambertian- 
scattering surfaces. The illuminated sliver of the planet at crescent phases represents a 
relatively small amount of surface area so that the illuminated disk at these phases can 
become dominated by large cloud features (or the glint spot), leading to large variability. 
Near full phase, the illuminated disk represents a large amount of surface area, averaging 
over clouded and non-clouded scenes, leading to relatively low variability. 

Our simulations without glint produce the same rise in variability into crescent phases 
as our simulations with glint, which indicates that variability at crescent phases is not com- 
pletely driven by glint. Thus, the trend of increasin g variability into cresce nt phases is not a 



clear indicator of surface oceans, as was proposed by lOakley fc Cash! (120091 ) . Any planet that 
can achieve sufficient contrast between bright and dark surfaces will produce a variability 
signal that increases into crescent phases, regardless of the presence of oceans. 

Our glinting model demonstrates a wavelength-dependent brightness excess over our 
non-glinting model, since some wavelengths are more sensitive to surface effects than others. 
The excess shrinks to less than 10% for the wavelength range 0.3-0.4 /im, where Rayleigh 
scattering obscures the surface. At wavelengths that correspond to relatively deep absorption 
bands, like the 1.4 fim water band, the excess shrinks to nearly zero because observations 
are mostly insensitive to the surface. 

At crescent phases, pathlengths through the atmosphere are relatively large and optical 
depths to Rayleigh scattering can be larger than unity even at longer wavelengths. This 
indicates that observations which aim to detect the brightness excess due to glint should be 
made at wavelengths in the near-infrared range. Earth's brightness drops by over an order 
of magnitude between 1-2 ^xm, arguing that searches for glint should occur below 2 /im for 
higher signal-to-noise ratio (SNR) detections. Since glint is a broad feature in wavelength 
space (it is the reflected solar spectrum, modulated by Rayleigh scattering, liquid water 
absorption at the surface, and atmospheric absorption), photometry can be used to detect 
glint provided that strong absorption features are avoided. 



5.2. Observing Requirements for Glint Detection 



The glint effect is strongest at NIR wavelengths, so that pairing JWST ( [Gardner et al. 
20061 ) with an external occulter (e.g.. lSoummer et al.ll2009l ; ICashll2006l ) would present a near- 
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term opportunity for the detection of oceans on extrasolar planets. Here we discuss optimal 
filter selection and inner working angle (IWA) and SNR requirements for glint detection at 
wavelengths accessible to JWST. 



The JWST Near Infrared Camera (NIRCam) (IHorner fc Riekdl2004l ) offers several medium 
and wide band filters suitable for glint detection, with the F115W, F150W, and F162M fil- 
ters (spanning 1.0-1.3 /im, 1.3-1.7 /im, and 1.55-1.70 /im, respectively) being most ideal. 
The F115W and F150W filters partially overlap water absorption features but experience a 
photon flux 3-4 times larger than the F162M filter, which does not span any strong water 
features and, thus, has an increased sensitivity to surface effects. 

Figure H] demonstrates a strategy that could be used to observe the brightness excess 
from glint. We show the 24-hour average brightness of Earth through the F115W NIRCam 
filter normalized to an observation at gibbous phase for several different orbital inclinations. 
Also shown is the SNR required to distinguish our glinting model from our non-glinting 
model at the l-a level assuming that high SNR observations (SNR=20) have been made 
at gibbous phase. For the inclination equals 90° and 75° cases, the glinting model shows a 
distinct leveling-off of its normalized brightness at crescent phases. Relative to its brightness 
at gibbous phase, the planet's brightness can remain roughly constant through certain phases 
due to the competing effects of the falling stellar illumination and the rising reflectivity due 
to glint. The glint effect becomes difficult to detect for the case where inclination equals 
60° because the planet never becomes a small enough crescent to produce a strong glint 
effect. Minimum SNRs between 5-10 are required for glint detection, depending on viewing 
geometry and telescope IWA. Note that roughly 25% of all planets discovered will have 
orbital inclinations between 75-105° and 50% will have inclinations between 60-120°. 

Observations of Earth-like exoplanets at crescent phases can be difficult due to IWA 
constraints. Vertical shaded regions in Fig. H] represent portions of the planet's orbit that 
cannot be observed for an Earth-twi n at different distances assuming an IWA of 85 milli- 



arcseconds (IBrown &: Soummerll2010l ). demonstrating that the glint effect could be detected 
for near edge-on orbits out to a distance of about 8 parsecs for this IWA. The measure of 
planet-star separation shown along the top of the sub-figures indicates that an IWA of about 
50 milli-arcseconds would allow for the easiest detection of the glint effect for planets within 
10 parsecs. Note that angular separation will scale inversely with distance to the system, 
and that the IWA constraints become significantly less strict for an Earth-like planet at the 
outer edge of the Habitable Zone. 
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6. Conclusion 

Our model successfully reproduces Earth's phase-dependent brightness from Earthshine 
observations. We have shown that glint increases Earth's brightness by as much as 100% at 
crescent phases, and that this effect is strongest in wavelength regions that are unaffected 
by Rayleigh scattering and atmospheric absorption. The glint effect may be detectable 
using JWST/NIRCam paired with an external occulter. Depending on viewing geometry, 
minimum SNRs of 5-10 are required for glint detection and an optimal IWA for detection is 
about 50 milli-arcseconds for an Earth-twin at 10 parsecs. 

This work was performed by the NASA Astrobiology Institute's Virtual Planetary Lab- 
oratory, supported under solicitation No. NNH05ZDA001C. 
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7. Tables and Figures 




Fig. 1. — A true-color image from our model (left) compared to a view of Earth from the Earth and Moon 
Viewer ( http: / /www.fourmilab.ch/cgi-bin/E arth /\ . A glint spot in the Indian Ocean can be clearly seen in 
the model image. 
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Fig. 2. — Simulation of Earth through a year. Black line corresponds to a model that includes glint 
while the grey line corresponds to a model that does not include glint. The left y-axis corresponds to 
the bell-shaped curves, demonstrating that Earth is brightest at full phase (orbital longitudes near 180°) 
and faintest near crescent phase (orbital longitudes near and 360°). The right y-axis corresponds to the 
bowl-shaped curves, where a perfect Lambert sphere would have a constant apparent albedo with phase. 
Variability at small time scales is due to Earth's rotation and time-varying cloud formations (noise is not 
included in simulations). Model "observations" are recorded every four hours, the system is viewed edge- 
on (i — 90°), and an orbital longitude o f 0° corresponds to January 1, 2008. Small stars are Earthshine 
measurements of Earth's apparent albedo (jPalle et al. 2003). Large circles are 24-hour average measurements 
of Earth's apparent alb edo recorded by the Deep Impact flyby spacecraft as part of NASA's EPOXI mission 



(jLivengood et al.l 12008). All data and model observations the wavelength range 0.4-0.7 /im. The bottom 



sub-panel demonstrates the brightness excess seen in the glinting model over the non-glinting model. 
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Fig. 3. — Variability in brightness for our glinting model (black) and our non-glinting model (grey), which 
are both shown in Fig. [2j Details are discussed in the text. 
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Fig. 4. — Earth's brightness through the JWST/NIRCam F115W filter (spanning 1.0-1.3 /jm) relative to 
its brightness at gibbous phase (135° and 225° orbital longitude) for orbital inclinations of 90° (top), 75° 
(middle), and 60° (bottom). Glinting model is in black and non- glinting model is in grey. Vertical shaded 
regions indicate the portions of the orbit for which a planet orbiting at 1 AU from its host star is withi n 85 
milli-arcseconds, which is a standard IWA for an occulter paired with JWST |Brown fc Soummerl[201f} ) . for 
a system at a distance of 5 parsecs (darkest), 7.5 parsecs (medium), and 10 parsecs (lightest). Planet-star 
separation at a distance of 10 parsecs is shown on the upper x-axis. The SNR required to distinguish the 
glinting model from the non-glinting model at the l-er level is shown along the right y-axis and corresponds 
to the dashed line. Observations have been averaged over 24-hour periods. 



